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Nuclear pairing interaction plays a crucial role in both macroscopic-microscopic and 
fully macroscopic descriptions of nuclei. In the present study we discuss different pairing 
interactions (monopole and S pairing forces) and the methods allowing for the particle 
number symmetry restoration in addition to the customary BCS treatment of pairing 
correlations in the context of o-decay half-lives for superheavy nuclei. The calculations 
are done in the macroscopic-microscopic framework for even-even nuclei with Z > 110. 

1. Introduction 

In the macroscopic-microscopic treatment of nuclear energy one needs to know the 
dependence of the total energy of the nucleus on deformation parameters (e.g., an 
elongation). In this method the energy is a sum of the macroscopic part (-E macro ) 
and microscopic energy consisting of shell (<5£shcii) and pairing (<5_E pa ; r ) energies 

V(def) = + C'-E'shell + SEpair . (1) 

The influence of the type of macroscopic energy models on both fission and 
on a-decay was discussed recently. 1 ' 2 For the macroscopic part we choose the 
Yukawa+ Exponential model (YpE). 3 The shell energy depends only on the single 
particle model - here Woods-Saxon potential. 4 The third component of the total 
energy is the pairing energy customary obtained with the monopole pairing force, 
i.e., the pairing force with constant matrix elements (g = const) in the BCS model. 
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As it is known, the usual BCS treatment of correlations leads to the particle num- 
ber non-conservation. To correct the solutions one may apply e.g., exact projection 
techniques, i.e., variation after projection method or use some approximative ap- 
proaches like solving Lipkin-Nogami equations. Both these methods are used here 
in the case of seniority pairing. 

In the following a short description of Lipkin-Nogami (LN) model and a variation 
after projection (VAP) method are presented together with the consequences of 
both these techniques on a-dccay energies (Q a -values) and a-half-lives (T a ) in the 
region of superheavy nuclei (Z > 110). 

2. Pairing models 

In this chapter we give a short review of the pairing models which will be used in our 
calculations. At each stage of the presentation we will also show some intermediate 
results concerning pairing properties of considered nuclei. 

State dependent S -pairing 

The 5 type pairing interaction induces the state dependence of the pairing gaps. In 
this case one has to solve a large number of BCS equations in order to obtain all 
the gaps in the pairing window (in this paper we solve A + 1 equations for neutrons 
and Z + 1 for protons). The strength of the 6- force which reproduces pairing gaps 
and masses in superheavy region of nuclei is 225 MeV fm 3 for both neutrons and 
protons. 2 In the case of monopole pairing we have used g parameters reported in 



Figures 1 and 2 show the Fermi level pairing gaps and pairing energies in both 
seniority BCS (BCSg) and state dependent BCS (BCS<5) models. The differences in 
pairing gaps are of the order of O.lMeV. The deviations in the pairing energy are 
in some cases larger than 0.5McV. 

Lipkin-Nogami method 

In the case of Lipkin-Nogami (LN) model one reduces the fluctuations of the par- 
ticle number by adding to the Hamiltonian H the terms — A A and the quadratic 
term — A 2 (A — (A)) 2 and one minimizes the average energy with respect to A. 
The procedure described in e.g., Ref. 6,7 gives the following expression for the new 
corrected energy 



Ref. 5 



E L n =E- A 2 (*|AA|*) 



(2) 



where 



(*|F(AA) 2 |*) 
(*|(AA) 2 |*) 



AA = A - (*|A|*) 



(3) 



N is the number operator and ^ is the BCS ground state. 
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Variation After Projection approach 



Variation after projection (VAP) or Full BCS (FBCS) is an exact particle number 
projection method 8,9 . Particle number projected wave function is 

* = C(fd( c" " 1 IJK + «vC44)*o , (4) 

where $o is the vacuum state for particles, hq the number of nucleonic pairs, and 
C the normalization constant 

|C| 2 = l/(-4^ 2 K, (5) 

where Rq is defined in Eq. (8). 

The energy of the system of N interacting fcrmions defined as 

E=mH\^)/m^}, (6) 

is given by 



\fc>o p,<?;p#' 



g pq u p v p u q v q Rl(p, q) ] . (7) 
The (hi, . . . , k n ) are the residues of contour integrals on the complex plane 

i r z n 

R%(k u ...,k N ) = —jdz— l J] (8) 

Here v k and u k are the parameters such that m 2 , + v\ = 1 and the number of pairs 
of nucleons reads 

n = ^M(k)/Rl (9) 

fc 

In order to calculate the residuum integrals we have used Ma and Rasmussen re- 
currence relations. 10 

The FBCS equations are highly nonlinear and hard to solve even in the case of 
seniority force gui = const. They read 

2{e k + A k )u k v k - A k (u 2 k - v%) = , (10) 



where 



e k = (e k -l/2 9kk )^-, (11) 
A k = 2^9kiuiv t — jtq — , (12) 



and 



Afc = m — E + l^ e ^i m 



2R> f—^ R' 

2>o 



1^ R%(k,l,m)- Rf(k,l,m) 

n / ,9lmUlVlU m Vm ■ 
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Fig. 3. o-decay energies for all models for isotopes of Z = 118. 



3. at half- lives 

Both the approximate LN and exact VAP projection methods have been applied 
to calculate macroscopic-microscopic masses, a-decay energies (Q-values) and half- 
lives (T a ) for even-even isotopes of superheavy nuclei: 112 < Z < 120. The a 
half- lives are determined from modified Viola-Seaborg formula. 11 In the following 
section we discuss the results of these calculations. 

The Qc-valucs 



Q a (Z, N) = [M(Z, N) - M(Z - 2, N - 2) - M(2, 2)]c 



(14) 



where M. is the nuclear mass, are depicted in Figure 3 for all considered models. 
The pattern of Q a vs. neutron number N is similar in all of the cases. However, 
in the vicinity of N — 162 and = 184 one observes some discrepancies which 
will influence the calculations of a-half-lives. This effect can be noticed on the 
Figure 4 where the logarithms of the half-lives (in years) are displayed. The largest 
differences occur in the vicinity of magic numbers where the standard BCSg and 
BCS<5 methods collapse yielding no superfluid solutions. 

Figure 5 shows the discrepancy between the half- lives (logarithms) vs. the 
neutron number A^ in the case of all of the models relative to the VAP results. 
The largest differences in logX^ are observed in the case of very heavy isotopes 
(N > 180). 
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Similarly, Figure 6 displays the correlations between the half-lives data for the 
log(T a /y) versus the data of the VAP model. The largest model discrepancies which 
reaches 4 orders of magnitude concern the long-living (log(T Q /y) > —8) isotopes- 
this corresponds to the region of superheavy nuclei which is of the main concern of 
both the experiment and the theory. 

4. Summary 

We have calculated a-decay half-lives using four pairing models: seniority (BCSg), 
(5-pairing (BCS<5), Lipkin-Nogami model and variation after projection method. As 
expected, both the kind of pairing interaction involved and the method of solution 
influence the a-decay half-lives. 

The largest discrepancies measured for 70 nuclei relative to the results of the 
VAP method are observed in the case of Lipkin-Nogami approach (an average 
deviation ctln = 0.717 MeV) and BCS with constant pairing interaction, where 
CBCSg = 0.654 MeV. Surprisingly, the (5-pairing BCS (BCS<5) model results are 
closer to the variation after projection(VAP) data in comparison to those obtained 
from the Lipkin-Nogami approach. In the latter case the mean deviation from VAP 
data a s = 0.438 MeV. 

The half-lives of long lived superheavy (N > 180) nuclei are determined with 
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the largest uncertainty. 

The problem of the description of the alpha decay of superheavy nuclei is still 
present in the realm of superheavy game and is a challenge for the theory to invent 
new, more powerful alpha decay models describing the phenomenon. 
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